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This paper is concerned with physics of the low energy singlet excitations found to exist below the 
spin gap in numerical studies of the Kagome lattice quantum Heisenberg antiferromagnet. Insight 
into the nature of these excitations is obtained by exploiting an approximate mapping to a fully 
frustrated transverse field Ising model on the dual dice lattice. This Ising model is shown to possess 
at least two phases - an ordered phase that also breaks translational symmetry with a large unit 
cell, and a paramagnetic phase. The former is argued to be a likely candidate for the ground state 
of the original Kagome magnet which thereby exhibits a specific pattern of dimer ordering with a 
large unit cell. Comparisons with available numerical results are made. 



I. INTRODUCTION 

Geometrically frustrated quantum magnets (GFQM) 
are a class of magnetic materials with various unusual 
properties at low temperatures. The combination of 
strong frustration and quantum effects makes them po- 
tentially a source of new strongly correlated physics. One 
common possible consequence of frustration is the macro- 
scopically large degeneracy of the classical ground-states, 
which makes them very sensitive to even weak ther- 
mal and quantum fluctuations. For example, the phe- 
nomenon of order-from-disorder may occur: the fluctua- 
tions may lift the initial degeneracy and yield a symmetry 
broken ground-state. Alternatively, the initial degener- 
acy may be completely lifted, in which case the result 
is a spin-liquid state. In the traditional point of view, 
the liquid state is a result of fluctuations strong enough 
to destroy the long-range order of the classical ground- 
state. However, in several GFQM, the starting classical 
ground-state is typically already disordered, so that the 
liquid state, eventually shaped by the fluctuations, may 
be of a different kind. In general, one can hope to find 
various exotic phenomena in these systems. 

The most studied GFQM are the Kagome and py- 
rochlore lattices of antiferromagnetically coupled nearest 
neighbor spins. They belong to a class of lattices which 
are composed of corner-sharing units, where every unit 
is a small frustrated spin system, like a triangle. Other 
examples include the triangular lattice frustrated by the 
higher order exchange processes, and the checkerboard 
lattice. 

This paper is concerned with the Kagome lattice spin- 
1/2 Heisenberg antiferromagnet (See Fig. QJ. Several 
groups have performed numerical calculations of spec- 
tra on finite clusters with up to 36 sitesA*^. The exact 
diagonalization2ii4 showed that the ground-state of the 
36-site Kagome lattice has short-ranged spin-spin and 
some other correlations. The spin-gap is finite, and the 
extrapolation to the thermodynamic limit suggests that 
it remains finite in the infinite lattice. It is, however, 
small, estimated around J/20, where J is the exchange 
coupling. Remarkably however, at energies below the 



spin gap there appear to be a large number of singlet ex- 
citations. The origin and nature of these states is unclear 
(see Refi^ for some interesting approaches). The num- 
ber of these states is found to be an exponential func- 
tion of the sample size: 1.15^ below the spin-gap, for 
even number of sites N, and there is no clear gap that 
separates them from the ground-state. There are some 
indications that these low-energy singlet excitations may 
also occur in other frustrated quantum paramagnets. In- 
deed, evidence for such excitations in the multiple-ring 
exchange triangular lattice spin-1/2 quantum antiferro- 
magnet has been reported in ReiX. This led to a sugges- 
tion that there was a new spin liquid phase characterized 
by having gapless singlet excitations, which was named 
the Type II spin liquid in Ref. 3 ^. 

There also are interesting experimental results on 
SrCrg p Gai2-9 P Oig (or SCGO, a bilayer spin-3/2 
Kagome magnet). Absence of long-range order at low 
temperatures^ accompanied however by a large entropy 
of low-lying excitations^ has been reported. Further- 
more, the heat capacity is virtually independent of mag- 
netic field, and is not thermally activated 9 -. These results 
suggest that the ground-state of SCGO is not magneti- 
cally ordered, and that the low-energy sector of excited 
states contains a (possibly gapless) band of singlet exci- 
tations, while spinful excitations are separated from the 
ground-state by a finite spin-gap. Despite the obvious 
differences between SCGO and the theoretical model of 
a single spin-1/2 Kagome layer, the similarities to the pic- 
ture obtained from exact diagonalization of the latter are 
striking. (Other unusual phenomena such as spin-glass 
like behavior in clean samples&iS and its speculated co- 
existence with a spin-liquid component 9 have also been 
suggested to happen in SCGO). 

These unusual properties of the Kagome Heisenberg 
magnet have motivated various theoretical efforts, but 
our understanding of the underlying physics is still very 
incomplete. Several controlled limits enable making some 
definite statements. One of them is the large- N SU(N) 
limit explored inii which gives highly degenerate mean- 
field solutions, and a spin-Peierls order as the 1 /N correc- 
tion. The ordered state maximizes the number of fluctu- 
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FIG. 1: (a) Kagome lattice; (b) dice lattice (dual to Kagome). 
The effective theory for the low-energy singlet degrees of free- 
dom will be described by a fully frustrated transverse-field 
quantum Ising model on the dice lattice. Frustration can be 
achieved in this model by giving opposite signs to the Ising 
spin interactions on the normal and thick bonds 



ating benzene-like hexagonal clusters, where three singlet 
bonds sitting on a hexagonal plaquette oscillate between 
two possible configurations. However, it remained un- 
known whether this ordered state would survive higher- 
order corrections. Another explored limit is the large-iV 
Sp(N) 12 . This generalization of the spin-1/2 problem 
retains the "spin-magnitude" as an independent parame- 
ter. For large spins the ground-state possesses long-range 
magnetic order (in a pattern dubbed the -\/3 x -\/3 struc- 
ture), and the elementary excitations are the spin waves 
associated with the broken spin rotational symmetry. For 
small spins, however, the ground-state does not break any 
symmetries, and the elementary spin carrying excitations 
are unconfined spin-1/2 bosonic spinons. 

An alternative, less controlled, but perhaps more phys- 
ical approach is to assume that the physics of the para- 
magnetic ground states of the Kagome magnet is de- 



scribed well by the short-ranged valence bonds, which can 
then be roughly approximated by the quantum dimers on 
the Kagome lattice. A quantum dimer model resulting 
from the overlap expansion of the valence-bond states 
was studied numerically by Zeng and Elseri Misguich 
and colaborators have considered two simple dimer mod- 
els: one at the Rokhsar-Kivelson points, and another 
offering more similarity to the overlap expansion (differ- 
ent resonant loops appear with different signs )i 4 . The 
former was exactly solved by mapping to the triangu- 
lar lattice Ising model, and both gave the dimer-liquid 
ground-states with large zero-temperature entropy. In 
a broader framework of frustrated quantum magnetism, 
Moessner, Sondhi, and various collaborators have stud- 
ied various related problems on different lattices, and dis- 
cussed connections between the dimer model and various 
gauge theories—. 

In this paper, we undertake a set of calculations that 
could potentially provide considerable insight into the 
low energy singlet excitations of the Kagome Heisen- 
berg antiferromagnet. It has been realized for some 
timet^!!^!^ that for a number of quantum antiferromag- 
ncts, the nature of the paramagnetic phases is closely 
related to (and essentially determined by) the properties 
of a much simpler model: the fully frustrated Ising model 
in a transverse field on a lattice that is dual to that of the 
original spin system. This connection (which we briefly 
review below) may be established in a number of differ- 
ent approaches. One is through a slave particle theory of 
such quantum paramagnets. The theory of fluctuations 
about the mean field state in such slave particle theo- 
ries is a gauge theory which for frustrated magnets^ is 
a Z2 gauge theory. This is then related by duality to the 
transverse field Ising model mentioned above. Another 
connection is through the quantum dimer description of 
the paramagnetic phases of the original spin system. As 
explained in Ref. 15 , in appropriate limits, the quantum 
dimer model may be related to the fully frustrated trans- 
verse field Ising model. 

The connection to the transverse field Ising model is 
particularly useful for issues related to the low energy 
singlet excitations of frustrated quantum paramagnets. 
First, by construction, the Ising model describes only the 
singlet sector of the original spin system, and is best to 
describe states below the spin gap. These are precisely 
the states one wishes to understand. Second, one can 
exploit the relative simplicity of the Ising model to per- 
form calculations that would be extremely difficult for 
the original spin system. For instance, the Ising model 
can be simulated by Monte Carlo without a sign problem, 
unlike the original frustrated spin model. This has been 
nicely demonstrated in the extensive work of Moessner, 
Sondhi, and coworkersi&iSi on such Ising models. 

It is therefore of extreme interest to study the particu- 
lar Ising model that one may guess describes the singlet 
sector of the Kagome Heisenberg antiferromagnet. This 
is a fully frustrated transverse-field Ising magnet on the 
dual lattice (known as the dice lattice; see Fig.^ 3 ). This 
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is an interesting model in itself, and appears not to have 
been examined before. As described in Appendix ]B\ a 
Landau-Ginzburg analysis (which has been useful^ in 
analyzing analogous frustrated quantum Ising magnets 
on other lattices) shows an infinitely degenerate set of 
zero energy modes in the Gaussian approximation. This 
renders the Landau-Ginzburg analysis less useful for this 
Ising model. Thus studying the fully frustrated quantum 
Ising model may also contain lessons of its own for the 
theory of frustrated quantum systems. 

The questions of interest are the following: Does this 
Ising model have a zero temperature ordered phase? 
We note that the combination of frustration and quan- 
tum effects could potentially completely destroy the or- 
dered phase even at zero temperature (as happens for 
instance for the transverse field Ising antiferromagnet on 
the Kagome lattic e 18 ' 19 ) . Generally, the ordered states of 
the Ising model describe, in the original Hciscnbcrg spin 
problem, "confined" paramagnets where the spin carry- 
ing excitations are (gapped) spin-1 magnons. The dis- 
ordered states of the Ising model describe fractionalized 
paramagnets where there are (gapped) spin- 1/2 spinous. 
Whatever the phase diagram of the dice Ising model, 
the nature of the excitations is an extremely appropriate 
question. Indeed, if this Ising model has a large number 
of low energy excitations, they would possibly correspond 
to the singlet excitations of the original Kagome lattice 
Heisenberg antiferromagnet. 

In this paper, we first argue that the fully frustrated 
transverse-field Ising model on the dice lattice supports 
(at least) two distinct phases. As usual, at large trans- 
verse fields f Section HV)l there is a gapped Ising param- 
agnet. However, its excitations are peculiar since they 
are either localized or extremely heavy. At small trans- 
verse fields (Section llllll there is also an ordered phase 
which breaks translational symmetry (albeit in a rather 
complicated pattern with a large unit cell). The ordering 
pattern is best described in the language of the quantum 
dimer model on the Kagome lattice to which the model is 
equivalent at small transverse fields: the dimers crystal- 
lize in a honeycomb structure of benzene-like resonating 
hexagons. The elementary excitations have a tendency to 
be localized (large effective mass), and they have a non- 
zero finite gap which nevertheless is very small compared 
to the natural energy scales of the model. We suggest 
that this phase is what is actually realized in the original 
Kagome Heisenberg magnet (see Fig.|2J), and not a gap- 
less spin-liquid. We argue that properties of this phase 
are consistent with results from the numerical studies, 
which in turn may suffer from serious finite-size effects. 
This then provides a qualitative picture of the singlet 
physics below the spin-gap observed by the exact diag- 
onalization. We also performed a simple Monte-Carlo 
simulation (Section [Vj) to show that the thermal fluctu- 
ations do not introduce order from disorder, and do not 
select the honeycomb structure out of many classically 
degenerate dimer coverings. In the appendices we pro- 
vide a discussion of various peripheral matters and details 




FIG. 2: Description of the suggested ground-state of the 
Kagome quantum antiferromagnet. Dimers formally repre- 
sent the frustrated bonds of the dual dice-lattice Ising model, 
but physically they represent singlet valence bonds. The 
translational symmetry is broken by a unit-cell with 36 sites. 
Dimers on the shaded hexagonal plaquettes resonate between 
the two possible configurations on the hexagon. Also, the 
six dimers around the central plaquette of every honeycomb 
cell resonate in a star shaped configuration. The lowest lying 
excitations are various excited states of the stars. 

of some calculations. 



II. MODELS 

We will try to understand the low-energy physics of 
the Heisenberg model on the Kagome lattice: 

H = jJ2Si-S j (1) 

(ij) 

where the sum runs over the nearest-neighbor sites. From 
the numerical exact diagonalization calculations we know 
that the ground state has no magnetic long range order. 
There is a non-zero gap for spin-carrying excitations. We 
will therefore focus on possible non-magnetically ordered 
phases of the model. A number of different arguments 
may be provided to relate the properties of such phases 
to those of the fully frustrated transverse field Ising model 
on the dice lattice. 

The quantum dimer model is a popular way to capture 
some essential physics of quantum paramagnetic ground 
states of spin systems. Moessner and Sondhi— have es- 
tablished a connection between the standard hard-core 
dimer model, and the fully frustrated Ising model in a 
(small) transverse field on the dual lattice. The connec- 
tion is made by noting that in zero transverse field, the 
ground states of the classical fully frustrated Ising model 
are (up to a global spin flip) in one-to-one correspon- 
dence with hard-core dimer coverings of the dual lattice. 
A small transverse field essentially introduces quantum 
resonances between dimer configurations, leading to the 
quantum dimer model. 
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A closely related earlier approach is through a slave 
particle mean field theory of the quantum paramagnetic 
phases of the Heisenberg model. As argued in Refill, for 
frustrated lattices the theory of fluctuations about the 
mean field theory is a Z2 gauge theory, where the gauge 
fields live on the links of the original (Kagome) lattice. 
For spin-1/2 systems, this gauge theory has non-trivial 
Berry phase terms (and was hence dubbed the "odd" Z2 
gauge theory in Ref>i£). Under a duality transformation, 
this odd Z2 gauge theory maps onto the fully frustrated 
Ising model in a transverse field on the dual (dice) lattice. 

Finally, as discussed ini£*2£, it is possible to directly 
formulate the Heisenberg spin model on any lattice as 
a theory of fermionic spinons coupled to an odd Z2 
gauge theory. In the spin-gapped quantum paramagnetic 
phases, the spinon fields may formally be integrated out, 
leaving a pure odd Z2 gauge theory as a description of 
physics below the spin gap. By the duality transforma- 
tion alluded to above, one again obtains the fully frus- 
trated transverse-field Ising model on the dual lattice. 
This approach gives a new flavor to it: the lattice does 
not have to be bipartite, and there is no need to rely on 
the dimer model from the beginning. We present a formal 
derivation and analysis in Appendix^ A rough formal 
estimate for the relative importance of various coupling 
constants is also provided. 

With these motivations, the starting point of this pa- 
per is the following frustrated Ising model in transverse 
field: 

H = -hJ2 zimvfv z m -X 3 £ vf 3 ~K 6 J2 < - 

(Irn) I3 ^6 

-K 3+3 £ «3--" (2) 

This Hamiltonian is defined on the dice lattice, dual 
to Kagome, and the frustration is realized through the 
ei m = ±1 factors, satisfying on each plaquette the fol- 
lowing relation: 

Y[e hn = 1 (3) 



The dice lattice, and a possible choice for e; m are shown 
in the Fig. \})p. The v x , v v , and v z operators are the 
Pauli matrices of the fluctuating Z2 vortex field which 
corresponds to the singlet degrees of freedom. From now 
on we will drop all terms denoted by the ellipses, and 
keep only the lowest order ones involving the isolated 3- 
coordinated sites ^3, 6-coordinated sites Iq, and the pairs 
of next-nearest-neighbor 3-coordinated sites (^37713). This 
theory is an effective theory of the Kagome Heisenberg 
magnet, describing the physics below the spin-gap. In the 
rest of the paper we study the properties of this model 
in different parameter regimes. As we formally argue in 
Appendix^ when obtained from the original Heisenberg 
model, it is natural to consider the limit of large h >> 
K3 6. This is also supported by comparison with the 
results of the numerical studies. Nevertheless we will 



study the model more generally and not just in the large- 
st limit. 



III. LARGE- /i LIMIT: VALENCE CRYSTAL 
PHASE 

Here we analyze the properties of the ground-state and 
excitation spectrum of the frustrated Ising model given 
by the Hamiltonian (J2J) in the h ^S> K3, Kq, -^3+3, . . . 
limit. We will build the perturbative expansion for the 
exact low-energy effective theory. 

The unperturbed Hamiltonian, with all K n equal to 
zero, is a classical Ising model on the fully frustrated 
dice lattice, and every spin configuration is an eigenstate. 
The ground-states have the minimal number of frustrated 
bonds (the bonds with positive energy: eimvfv^ = — 1). 
In trying to construct them, we have to leave at least 
one bond frustrated on every plaquette: otherwise, the 
product of four bond energies around a plaquette could 
not be negative, as required by the frustration condition 
©. It is convenient to return to a dimer representation 
of these states on the original Kagome lattice: we put 
a dimer on every Kagome bond which intersects a frus- 
trated dice bond. This representation is accurate up to 
global spin flip of the Ising spins. There will be an odd 
number of dimers emanating from every Kagome site, 
and the ground-states will be given by the various hard- 
core dimer coverings. 

The perturbations K n introduce quantum fluctuations 
between different dimer coverings, and lift degeneracy of 
the unperturbed ground-states. The effect of a vf opera- 
tor in the dimer picture is to toggle the dimer occupancy 
on all bonds of the Kagome plaquette which corresponds 
to the dice site I. Since we want to construct the ef- 
fective theory, we need to find the combinations of vf 
operators which transform one hard-core dimer covering 
into another, possibly through some virtual states. Such 
processes can be described as dimer motion along the 
flippable loops, and we will call them the loop flips (see 
Fig- El for explanation). It is convenient to introduce the 
following terminology: if a plaquette carries no dimers 
on its bonds, we will call it a defect plaquette, and if it 
carries a flippable loop on its bonds, we will call it a per- 
fect plaquette. A plaquette can be perfect only if it has 
an even number of bonds. Of special interest are the ele- 
mentary flippable loops which enclose only one hexagonal 
plaquette, and they are explained in Table [I] 

The lowest order terms in the perturbation theory are 
the single plaquette flips. The K3 and ^3+3 terms cannot 
connect between two hard-core dimer coverings, so that 
to the lowest order the effective theory is: 

H = -Nh - K 6 W 6 (4) 

where N is the number of Kagome sites, and Wq is the 
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(a) (b) 

FIG. 3: (a) An example of a flippable loop. The flippable 
loops are tangential to the dimers, so that they cannot contain 
only one end of any dimer. (b) An example of a loop-flip 
between the two possible hard-core dimer arrangements on it. 
This loop contains a hexagon and a bowtie pair of triangles. 
The loop is flipped through a virtual state by the successive 
operation of — K^vf and — K3+3V2V3 , where 1,2,3 refer to the 
dice lattice sites which correspond to the depicted Kagome 
plaquettes. 




(a) (b) 

FIG. 4: Two perfect hexagons cannot be nearest neighbors, 
and the perfect hexagon flip does not affect length of any other 
elementary flippable loop. Consider a perfect hexagon 2 being 
flipped from the configuration (a) to (b). The flippable loop 
on the hexagon 1 has to go through the sites A, B, and C, 
and therefore pass through the shaded bonds in order to be 
tangential to the dimers. As a consequence, it has to include 
at least one triangle, so that it cannot be a perfect hexagon, 
and also, its length is not affected by the flip on the hexagon 
2. 



kinetic energy operator of the perfect hexagons: 



W 6 



E 

o 



+ 



(5) 



The sum runs over all hexagonal plaquettes of the 
Kagome lattice. The quantum fluctuations created by 
the kinetic term yield the ground states in which the 
number of resonating perfect hexagons is maximized. 
This is an exact statement, since all perfect hexagon 
flips commute with one another: as illustrated in Fig. 0] 
the perfect hexagons can never be nearest neighbors, 
and therefore cannot affect one another. In order to 
gain more insight about these states, we want to re- 
call some observations from Referenced andS. The 







TABLE I: The elementary flippable loops: (a) flippable 
hexagon, (b) flippable 8-bond loops (rhombus, arrow and 
trapeze, from left to right), (c) the star, (d) flippable 10-bond 
loops. The elementary flippable loops by definition enclose 
only one hexagon. For arbitrary dimer covering, a unique ele- 
mentary flippable loop can be found on every hosting hexag- 
onal plaquette: it goes through all the hexagon sites, and 
includes those surrounding triangles which hold a dimer on a 
bond that does not belong to the hexagon. The hexagon in 
(a) is perfect, while the hexagon in (c) is a defect hexagon. 
The length of the elementary loop is directly related to the 
number of dimers sitting on the hexagon, which is equal to 
the number of defect triangles around the hexagon: 3 in (a), 
2 in (b), 1 in (d), and none in (c). 



total number of dimers in a hard-core dimer covering 
on the Kagome lattice is N2 — N/2, and the num- 
ber of triangular plaquettes is N A = 2N/3, so that 
N2 = 3N&/4 = (1 — l/4)A r A- Since a triangle can carry 
at most one dimer, we see that one quarter of all tri- 
angles are the defects in any hard-core dimer covering: 
N/\d = N&/4. Next, we note that every perfect hexagon 
has exactly three neighboring defect triangles, and since 
no two perfect hexagons can be neighbors, those defect 
triangles are not shared between them. It follows from 
this that the total number of perfect hexagons has an up- 
per bound: N 0p s£ N Ad /3 = N A /12 = N /6. The max- 
imum possible density of perfect hexagons is one per six 
hexagonal plaquettes, and it can be achieved in a variety 
of ways. In Fig. we show two characteristic possibili- 
ties, (a) the honeycomb, and (b) the stripe state. In gen- 
eral, these states are constructed by placing the perfect 
hexagons as close as possible to each other. The closest 
they can be is the next nearest neighbors, provided that 
between them is another hexagonal plaquette, and not a 
bowtie pair of triangles (x), because in the latter case 
there would have been a site (the center of the bowtie) 
involved in no dimers. This rule allows one to arrange 
perfect hexagons in strings which may be straight, bent 
at an angle of 120°, or forked into two new strings at the 
120° angles. The stripe state is an example of straight 
strings, while the honeycomb state has strings forking at 
each perfect hexagon. 

Before we proceed with the next order of the pertur- 
bation theory, we need to make some additional remarks. 
If we look at the elementary flippable loops realized on 
the various hexagonal plaquettes in Fig. we observe 
that between every two closest perfect hexagons there is 
an 8-bond flippable loop, right on the sides of the strings 
there are only 10-bond flippable loops, and in the case of 
the honeycomb state, there is a 12-bond star-shaped flip- 
pable loop sitting at the center of every honeycomb cell. 
Also, one never finds an arrow-shaped 8-bond flippable 
loop between two perfect hexagons. These are quite gen- 




(a) 




(b) 



FIG. 5: Two hard-core dimer patterns that maximize the 
number of perfect hexagons: (a) honeycomb pattern; (b) 
stripe pattern. The perfect hexagons are shaded to guide 
the eye. Note that the 8-bond fiippable loops appear only 
as "connections" between the perfect hexagons, and the 10- 
bond fiippable loops touch exactly one perfect hexagon. The 
honeycomb pattern has the 12-bond fiippable loops, the stars: 
they sit inside the honeycomb cells. 



eral features of the states with the maximum number of 
perfect hexagons, which we explain in more detail in the 
Fig. El 

At the second order of perturbation theory, we need 
to include the combinations of two K$, Kq, and -K3+3 
flips. One of them is the flip on the arrow shaped 8-bond 
fiippable loop (a hexagon, and a bowtie: see Fig[3t>): 



K 6 K 3+3 w (ar.) 



2h 



W> 



(ar.) 



E 

o 



+ 



(6) 
rotat. I (7) 



A B 

C D 

(a) 




All other allowed combinations flip the same plaquette 
twice, leaving the dimer configuration unchanged. A 



(b) 



FIG. 6: 

(a) The hexagonal plaquette between a pair of perfect 
hexagons always hosts an 8-bond fiippable loop. Once the 
perfect hexagons 1 and 3 are placed, the only way for all of 
the sites A,B,C, and D to be involved in dimers is to pair 
A,B and C,D. Then, the hexagonal plaquette 2 carries two 
dimers on its bonds, and therefore hosts an 8-bond fiippable 
loop. However, it cannot be an arrow-shaped loop (see Ta- 
ble EJ , because both shaded bowties must have one normal 
and one defect triangle due to the perfect hexagons, and an 
arrow-shaped loop must contain a bowtie with two normal 
triangles. 

(b) Only 10-bond fiippable loops surround the perfect 
hexagon strings. Consider a string, and a neighboring 
hexagon 1. One cannot put there a 6-bond fiippable loop, 
because this would be a perfect hexagon next to another per- 
fect hexagon from the string. A star-shaped 12-bond fiippable 
loop is not an option either, because this hexagon always has 
one neighboring defect triangle from the string. There is only 
one way to place two dimers on the bonds of hexagon 1 in at- 
tempt to create an 8-bond fiippable loop on it. However, this 
requires the site A to pair with either B or C, which in turn 
makes it impossible for another string of perfect hexagons to 
take its normal place. At least a spot is lost for a perfect 
hexagon, and this costs energy K$\ Similar is true for the 
corners of the outward bending strings, while this issue does 
not occur for the inward bending, as in the case of the hon- 
eycomb state, since 10-bond loops are the only choice there. 

Relative positions of various kinds of elementary fiippable 
loops can also be established by using the relation between 
their length and the number of defect triangles around them 
(Table [IJ • When the number of perfect hexagons is maxi- 
mized, all defects are placed around them (NAd = 3Nq p ), 
and the length of any fiippable loop depends on the number 
of its perfect hexagon neighbors. 
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hexagonal plaquette can be flipped twice if it is not a 
flippable loop: it has to be the part of an elementary 
fiippable loop with 8, 10, or 12 bonds. As the potential 
energy associated with the flip is different in each case, 
we need to know the numbers U n of the n-bond flippable 
loops in the dimer configuration in order to write the 
appropriate contribution in the effective theory: 



AH, 



(i) _ 

Ox2 ~~ 



K 2 

6 TT 



K 2 
6 TT 

8h Ul 



K 2 

6 TT 

Y2h Ul2 



(8) 



The number operators U n are explicitly written in the 
Table HH 

We treat the triangular plaquette double flips similarly. 
Every triangle can be cither a defect, or normal in any 
dimer covering. The defects contain no dimers, while the 
normal triangles contain one dimer. If their numbers are 
N Ad and N A respectively, then the contribution of the 
triangle double flips is: 



AH 



(i) 

Ax2 



K 2 
~-N 
6h 



K 2 
2h 



(N A -N Ad ) 



(9) 

h 18 ( ' 



We have used the identities N Ad = N A /A, and N A = 
2N/3, where N is the number of Kagome sites. Finally, 
it remains to consider the double flips of bowties (cxi). 
They can contain either one defect triangle or none. For 
every defect there are three bowties containing it, while 
the total number of bowties in the lattice is equal to the 
number of Kagome sites N. Therefore, the contribution 
of the bowties is: 
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We see that the triangle and bowtie contributions reduce 
to mere constants for arbitrary dimer covering. Then, 
it pays to investigate closer the contribution of double 
hexagon flips. First, we want to express the number of 
defect triangles in terms of the U n operators. This can 
be achieved by noting that every type of elementary flip- 
pable loop has a fixed number of both kinds of triangles 
(see TableHJ: N Ad = (3U 6 + 2U S + U 10 )/3, where the fac- 
tor of 1/3 corrects the over-counting of triangles shared 
between the elementary flippable loops. Using this, we 
can eliminate U$ and [/ 10 from (JHJ): 
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Now, we can add (0 E9 EH i an d write the effective 
theory to the second order: 
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The ground-state of this Hamiltonian can be obtained by 
the second level perturbation theory in which the unper- 
turbed Hamiltonian is given by g|. The W 8 (or ° term 



slightly spreads the ground-state wavefunction from the 
sharp honeycomb or a stripe-like state, but the correc- 
tion to the ground-state energy due to this term ap- 
pears only at the higher orders: both the honeycomb and 
stripe states cannot contain the arrow-shaped 8-bond 
flippable loops, and even if they could, the flip would 
destroy the two neighboring perfect hexagons, so that 
(0|W^ 8 (ar ) |0) = 0. Therefore, the energy shift is domi- 
nated by the potential energy part of 112fl involving U n 
operators. This part fully commutes with ijlfl: the per- 
fect hexagon flips Wq cannot change the length of any el- 
ementary flippable loop (see Fig. QJ , and the U operators 
simply count the number of flippable loops with the given 
length. Only the double hexagon flips Hill select the ac- 
tual ground state: when the number of perfect hexagons 
is maximized, the Uq operator behaves as a number, tak- 
ing the value Nq/6 = N/18, so that the potential energy 
is controlled by the number of the 12-bond star-shaped 
flippable loops. The honeycomb state (Fig. |Sp,) maxi- 
mizes their number, and therefore the quantum fluctua- 
tions select it as the ground-state. The degeneracy of this 
ground-states is exponentially large in the system size, at 
this order of the perturbation theory. It comes from the 
freedom to flip any star-loops without energy cost. 

Various new kinetic and potential terms appear in the 
effective dimer model at the higher orders. They further 
spread the ground-state wavefunction, but its main com- 
ponent remains the honeycomb structure. The ground- 
state degeneracy is finally lifted when the star kinetic 
terms appear (only the star-flips do not destroy the per- 
fect hexagons): 
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(13) 



In terms of K n /h, this first happens at the fourth order 
by combining the flips on one hexagon and three bowties: 
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However, since in principle it might happen that -K3+3 <C 
K 2 /h, the dominant star kinetic term may occur when 
a hexagon flip is combined with six triangle flips: 
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In any case, the ground-state will be the honeycomb state 
with resonating perfect hexagons and stars, while the 
lowest lying excitations will be gapped. The gap is very 
small, of the order of K 6 K% +3 /h 3 or K e K% / h 6 . The num- 
ber of such excitations is exponential in the system size. 
The only remaining degeneracy is 12-fold, due to the bro- 
ken translational symmetry (the honeycomb unit-cell has 
12 hexagonal plaquettes). 

It is also interesting to understand this ground-state in 
the dual Ising picture. As explained at the beginning of 
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TABLE II: The flippable loop number operators: E7 n counts 
They correspond to the potential energy in the dimer model. 

this section, a dimer on the Kagome lattice represents a 
frustrated bond of the dual Ising model on the dice lat- 
tice. Every Kagome dimer covering corresponds to two 
Ising spin configurations on the dice lattice, related to 
each other by the global spin-flip. If the Kagome dimers 
are flipped along a flippable loop, the corresponding ef- 
fect in the dual picture is the simultaneous flip of all dice 
Ising spins which sit "inside" the Kagome plaquettes en- 
closed by that flippable loop. Now we can translate the 
description of the ground-state. Every resonating perfect 
hexagon is a flippable spin on a 6-coordinated site of the 
dice lattice, in the state of equal superposition of "up" 
and "down". Every resonating star-shaped flippable loop 
is a flippable cluster of seven Ising spins coherently fluctu- 
ating between two states of defined spin orientation (one 
spin is on a 6-coordinated site, and the other six are on 
the surrounding 3-coordinated sites). However, certain 
dimers are static in the ground-state: two of them reside 
between every pair of neighboring perfect hexagons. In 
order to describe them in the dual language, it is suffi- 
cient to arrange the corresponding Ising spins on the dice 
lattice in some appropriate static configuration. There- 
fore, the ground-state breaks the global spin-flip symme- 
try of the dual Ising model, with 3/4 of all spins assum- 
ing a fixed orientation, and 1/4 of spins fluctuating. The 
translational symmetry is broken only by the arrange- 
ment of frustrated bonds, and locations of the fluctuating 
spins; there need not be actual long-range order in terms 
of the orientation of non-fluctuating spins, since this de- 
pends on a relatively arbitrary assignment of ei m = ±1 
in (J2J. In the dual language, this state is formed first by 
minimizing the number of frustrated bonds, and then by 
maximizing the number of flippable spins which gain the 
kinetic energy from the transverse field. 

Since this is a valence-bond crystal phase which breaks 
the spin-flip symmetry of the dual Ising model, it is also 
a confined phase. The elementary spinful excitations are 
gapped spin-1 magnons. We also note that this phase 
is stable against the weak fluctuations of the large loops 
that we have ignored from the Hamiltonian @. The 
stability has been demonstrated in the perturbation the- 
ory. 



number of n-bond flippable loops in the dimer configuration. 



At the end, we ask what changes if the system be- 
comes finite. As we have seen, the infinite system breaks 
translational symmetry. The ground-states are 12-fold 
degenerate, and related to one another by translations. 
Roughly speaking, we can even regard the lowest lying 
excited states as 12-fold degenerate. If the lattice be- 
comes finite, this remaining degeneracy will be lifted, and 
no symmetry will be spontaneously broken. Let us apply 
the perturbation theory developed in this Section to a fi- 
nite sample with TV sites. The 12-fold degenerate ground 
(and excited) states, shaped at the low orders, can be 
eventually mixed only by high order perturbations, which 
are essentially just translations. In the dimer language, 
they correspond to dimer flips on large flippable loops 
which enclose a finite fraction of all plaquettes of the sam- 
ple. Therefore, their order scales as N, and the energy 
scale associated with them is exponentially small. This 
is enough to completely mix the states with the same 
energy, but there will be no mixtures of the states with 
different energies. As a consequence, the spectrum of the 
finite system will be very similar to that of the infinite 
system, and just slightly more dispersed. For example, 
let us consider a sample with N — 36 sites. It takes en- 
closing at least ~ 16 plaquettes by the translational flip- 
pable loops, so that the translational degeneracy is lifted 
at ~ 16 th order. The exact count of the lowest lying exci- 
tations depends on the boundary conditions, but we can 
still demonstrate that this number is not small. If the 
boundary conditions allow the honeycomb pattern, then 
there will be at most one star and two perfect hexagons 
for N — 36. Some basic excited states can be obtained 
by exciting the star and two perfect hexagons in various 
combinations (2 3 total, including the ground-state), and 
each is 12-fold "degenerate" in the thermodynamic limit. 
The other can be obtained by departing from the honey- 
comb pattern, which may require removing one or both 
perfect hexagons. The "degeneracy" is lifted in this fi- 
nite system, and these states form bands. Their number 
is greater than 2 3 x 12 = 96. This is not inconsistent 
with the number obtained from numerics: 1.15 ~ 153 
below the spin-gap. We also note that N = 36 is small 
enough to miss the spinon confinement. Confinement is 
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noticeable only at distances much larger than size of the 
honeycomb pattern unit-cell, so that seeing it would re- 
quire at least N = 72 sites. 



IV. SMALL-ft LIMIT: FRACTION ALIZED 
SPIN-LIQUID PHASE 

The limit h -C K 3 , -K3+3, Kq. . . is convenient to ana- 
lyze directly in the frustrated Ising model. If h vanishes, 
then all Ising spins in (J2J are aligned with the transverse 
fields in the x-direction, making the ground-state com- 
pletely disordered and uncorrelated in terms of vf. An 
elementary excitation is created when one spin is flipped 
against the transverse field. These excitations are Z 2 
vortices, or visons according to (|A7I) . and they are local- 
ized and gapped. For finite h, the visons can in principle 
hop between sites and lower their gap by acquiring ki- 
netic energy. However, for small h the gap is guarantied 
to persist, and the ground-state remains disordered, and 
unique. 

In the following we will assume that K 3 and K§ are 
positive, and we will ignore ^"3+3 and other terms for 
simplicity. One way to study the properties of the excited 
states is to consider perturbatively the effective theory 
for one vison. At the lowest order, this effective theory 
is simply the nearest neighbor hopping Hamiltonian: 

H = -hJ2 (\l)eim(m\ + h.c^j + 

(lm) 

+2K 3 J2\ l 3)(h\+2K 6 J2\le)(k\ (16) 

h h 

where ^3 and Iq are the 3 and 6-coordinated sites of the 
dice lattice, and |^) = vf\Q). It can be easily diagonal- 
ized in the momentum space by working with the 6-site 
elementary cells of the tile in Fig.Q] It was shown in the 
reference^ that for K 3 = Kq the energy spectrum of this 
model is completely dispersionless and divided into three 
macroscopically degenerate levels. Remarkably, this re- 
mains true for arbitrary values of K 3 and Kq : 

E x = 2K e 

E 2 = 2K 3 - V6h (17) 
E 3 = 2K 3 + V6h 

Therefore, the visons are localized at the lowest order of 
perturbation theory, in the similar way to a single elec- 
tron in the magnetic field (this analogy does not generally 
hold). A natural question to ask is whether this localiza- 
tion persists to higher, or maybe all orders of perturba- 
tion theory. We try to find an answer by considering the 
time-ordered Green's function in the interaction picture: 

ir(! m . f t s (Q|r^)?4(US|o) (m 

iLr(l, m\ t h t m ) (01510) 




(a) (b) (c) (d) 

FIG. 7: Examples of the paths between the sites P and Q: (a) 
simple path; (b)-(d) complex paths. The thin line-doublets 
around the bonds in (c) and (d) mean that there are loops 
going through those bonds twice. 

where T is the time-ordering operator, and: 
S = Texp(-i [ &tH'{t)\ 

J — GO 



v z (t) = e lHot v z (0)e- lHot (19) 

The perturbation H' is the vison-hopping part of J3J), 
while Hq is the vison potential energy. The expansion of 
S in (|18f) generates many vacuum expectation values of 
the products of v z operators, which are then integrated 
over the internal time variables. These vacuum expecta- 
tions are easy to calculate for any given set of time mo- 
ments: the v z (t) operators can be evolved to t = when 
they simply flip a spin, while the evolution parts e °* 
only introduce a phase factor, since they always act on 
a Hq eigenstate. We immediately see that a vacuum ex- 
pectation which appears in l|18|l can be non-zero only if 
the product of v z is made from a path of bonds bridging 
between the sites I and to. We only need to turn our 
attention to the connected clusters of bonds. 

Consider a particular pair of sites P and Q, and 
the n th order term in the perturbative expansion of 
iG(P,Q;tp,tQ). Its value, proportional to h n , is con- 
tributed separately by all paths of the length n between 
P and Q. The two paths cancel out each other if they 
enclose an odd number of rhombic plaquettes: the to- 
tal "flux" through the loop formed by them is tt, that 
Ilioop e i"i = — 1- Let us distinguish the two types of 
paths: simple and complex. The simple paths never visit 
any site more than once, while the complex paths visit at 
least one site more than once. The complex paths can be 
obtained from the simple ones by adding the loop seg- 
ments to them, where a loop may take some bonds an 
even number of times. We illustrate this in the Fig. [7] 

For any given simple path between some two distant 
sites P and Q, we can construct a path that cancels it 
out. In demonstrating this, we also explore what the 
sufficient distance between the P and Q is. (a) Con- 
sider first the simple paths which contain only the three 
6-coordinated sites depicted in the Fig. |SJi. If the path 
goes through the sites A\ — B\ — A2, and never visits 
the site B2, then it can be canceled by the path which 
goes through A\ — B2 — A2, and continues beyond as the 
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(a) 



(b) 



(c) 



FIG. 8: If the paths contain only the marked 6-coordinated 
sites, they would have to contain also all of the marked 3- 
coordinated sites in order not to be cancelable (see the text). 



original one. The only way to avoid the path cancella- 
tion is to visit all the depicted B-sites, and this can be 
done only by choosing the path's end points among them 
(otherwise, more 6-coordinated sites would belong to the 
path) . We see that in this case the path's end-points are 
the next-nearest- neighbors, (b) Next, consider the paths 
which contain the four 6-coordinated sites depicted in 
the Fig. |SJd. For the same reason as before, all the B- 
sites must be included in the path, otherwise there will 
be another path which cancels it out. But this time, in 
attempting to do so we would have to introduce another 
6-coordinated site into the path (the remaining neigh- 
bors of the depicted B-sites - analyzed as the case c), 
since the path-ends can take only two out of three B- 
sites left "outside" of the straight segment of the path. 
We can already see that all the simple paths with more 
than three 6-coordinated sites arranged in a chain can- 
cel out. This chain may, of course, bend, but somewhat 
special situation occurs when a triplet of 6-coordinated 
sites (Ai, A 2 , A 3 ) of the path sits on the three touch- 
ing plaquettes, as shown in the Fig. QJ. (c) Again, all 
the B-sites must be visited, or the path is canceled out, 
and one of them must be a path's end-point if no other 
6-coordinated sites may belong to the path. Since all 
B-sites are now the next-nearest-neighbors, we may try 
to choose another 3-coordinated site further away as the 
path's end-point, and construct the path like the one 
plotted with the solid line in the Fig-Efc- But, again there 
is a path which cancels it out, and it is plotted with the 
dashed line. By trying to add more 6-coordinated sites 
to the case c) we would only have too many B-sites to 
cover as before. In conclusion, all the simple paths whose 
end-points are separated more than the next-nearest- 
neighbors must cancel out. 

Many complex paths are also canceled - namely, when- 
ever we can construct another path of the same topology, 
and between the same end-points, which encloses an odd 
number of plaquettes with the original one. However, for 
some paths this construction is not possible, and a differ- 
ent approach is needed in order to see whether they are 
canceled out or not. An example of such a path which 
appears at the 12 th order of the perturbation theory is 
shown in the Fig. [7|i. The complex paths between the 
sites further away than the next-nearest-neighbors ap- 
pear only above the 6 th order, so that below this order 



the Green's function between distant sites (and arbitrary 
times) strictly vanishes, and the visons appear strictly 
localized. 

At present, we cannot tell whether the vison localiza- 
tion persists to all orders. If not, they would have an ex- 
tremely large inertia in the small- h limit, and this would 
be a novel mechanism for creation of the large quasi- 
particle effective mass. The possibility of vison localiza- 
tion is also hinted by the Monte Carlo calculations pre- 
sented in the Section ^ Namely, the spatial vison- vison 
correlations effectively vanish at separations beyond the 
next-nearest neighbors. We would like to contrast this 
with another case of a fractionalized spin-liquid on the 
Kagome lattice: in a particular easy-axis limit explored 
in Ref 2^ the vison correlations are short-ranged, but they 
exhibit a clear exponential decay over a significant range 
of separations. 

Finally, we note that the phase obtained here is also 
the unconfined phase. The elementary singlet excitations 
are gapped visons, while the elementary spin-carrying 
excitations are gapped S = 1/2 spinons. 



V. MONTE-CARLO RESULTS 

In order to gain some more insight into the frustrated 
transverse field Ising model, we attempted to study it nu- 
merically using the classical Monte Carlo technique. The 
quantum Ising model (J2J translates to the 3-dimcnsional 
classical Ising model on the time-stacked dice-lattice lay- 
ers: the corresponding Boltzmann weights are all positive 
and hence the model can be simulated by Monte Carlo 
without a sign problem. We reliably measured the spin- 
spin correlations in the disordered phase, and answered 
the question of whether the thermal fluctuations alone 
can yield order-from-disorder like the quantum fluctua- 
tions. However, we could not truly explore the complete 
phase diagram and the nature of the phase transition 
with this simple classical technique; for those purposes, 
a more sophisticated Monte Carlo technique is needed. 

The model we study is given by the action analogous 
to 10: 



T lm 

-K 3 ^2v htT v h . T+1 - K 6 ^2vi 6 , T vi 6 , T+1 (20) 

We have to regard this problem as a classical one, be- 
cause the quantum limit h — > 0, and K^^ — > oo is not 

accessible. If K 3 = K 6 — K, two regimes emerge, de- 
picted in the Fig. [5] In order to obtain this diagram, we 
scan along K oc h lines, and measure the heat-capacity, 
using h as the inverse "temperature" /3; then we record 
position of the heat-capacity peak, and join the peak po- 
sitions [h, K) from different scans into a crossover line 
between the two regimes. 
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For small h and K a disordered regime occurs. It 
is characterized by seemingly vanishing spatial spin-spin 
correlations beyond the nearest-neighbor sites, as shown 
in Fig. I1UI while the time correlations decay exponen- 
tially. This regime corresponds to the disordered phase 
and the small- h limit of (J2J). The numerical result for 
the spatial correlations is consistent with the statement 
made in the Section Hvl that the visons may be localized 
in the vicinity of the next-nearest-neighbor sites. 

For large h and K, a "dimerized" regime occurs: if 
the frustrated dice bonds (with positive bond energy) 
are plotted as dimers on the Kagome lattice (wherever a 
Kagome bond intersects a frustrated dice bond), then 
a hard-core dimer covering is obtained. Our present 
method was not able to establish whether the crossover 
line between these two regimes is actually a phase transi- 
tion line or not. However, based on our analytical study 
we expect that this is the phase transition. Since large 
coupling constants yield the dimerized regime in which 
all dimer coverings have the same action, the remaining 
phase properties can be decided only by entropy, which 
is a manifestation of quantum fluctuations. In the Ising 
language, once the action is fixed to its lowest value, the 
maximum entropy state has the largest number of flip- 
pable clusters of spins (free to fluctuate without chang- 
ing the action). In the dimer language, the honeycomb 
pattern will be selected, in a mechanism essentially anal- 
ogous to that obtained in the Section ITTT1 Indeed, every 
perfect hexagon of dimers represents a flippable spin on 
the dice-lattice. Once their number is maximized by the 
entropy, the only other elementary flippable loops which 
could possibly resonate without action cost are the stars, 
and to maximize their number, the perfect hexagons must 
arrange in the honeycomb pattern. Unfortunately, our 
simulation could not reach equilibrium in this regime, at 
least due to very slow dynamics at large coupling con- 
stants. 

Since the quantum fluctuations are capable of produc- 
ing the order-from-disorder, it is natural to ask whether 
the thermal fluctuations are capable too. If we set the 
transverse field couplings in J2J to zero, we obtain the 
classical two-dimensional Ising model on the fully frus- 
trated dice lattice. In this case, the Monte Carlo eas- 
ily reveals that there are no phase transitions, measured 
through the "heat capacity" , and the model is always in 
its disordered phase. Therefore, the thermal fluctuations 
alone cannot introduce a long-range order, even when h 
becomes large and the "dimerized" regime occurs. The 
presence of the transverse fields is crucial for the entrop- 
ical selection of a long-range ordered state in the Monte 
Carlo. 



VI. DISCUSSION: IMPLICATIONS FOR THE 
KAGOME QUANTUM HEISENBERG MAGNET 

The analysis in the preceding sections has established 
some properties of the fully frustrated transverse field 
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FIG. 9: The regimes of the (2+l)D classical Ising model on 
the fully frustrated dice lattice in I2UI . K = K3 = K§. The 
crossover line was obtained by joining the points (h, K) at 
which the heat-capacity is peaked for several fixed ratios K/h. 
In defining the heat-capacity h was treated as the inverse 
temperature j3. 



Ising model on the dice lattice. Specifically, we analyzed 
the two limits of small and large transverse fields K n in 
(J2J) , and described the two distinct phases that obtain in 
those limits. In this Section, we consider the lessons for 
the original Kagome Heisenberg magnet. The disordered 
phase of the Ising model (h ^ K n ) has a unique ground- 
state and a clear gap in its spectrum. This will describe 
spin liquid phases of the original Kagome magnet. On 
the other hand, the honeycomb dimer phase (h K n ) 
breaks translational symmetry, but with a very large unit 
cell. This phase has a number of interesting properties. 
First it clearly describes a valence bond crystal phase of 
the original Kagome magnet with a large unit cell. The 
spinful excitations will be gapped spin-1 magnons. How- 
ever, below the spin gap this state has a large number 
of low-energy excitations with a very small but non-zero 
gap. The number of low-lying excitations that we can 
account for is exponential in the system size. However, 
their exact number below the spin-gap is hard to esti- 
mate in our approach. The numerical calculations^ on 
the Kagome Heisenberg magnet revealed a large num- 
ber of seemingly gapless singlet excitations in the 36-site 
sample of the Kagome quantum magnet, and a finite- 
size scaling analysis suggested that in the thermodynamic 
limit the spectrum could be gapless. We suggest that the 
actual ground state in the thermodynamic limit displays 
the honeycomb dimer ordering present in the dual Ising 
model at small transverse fields. At the end of Section 
IIIII we argued that our suggestion is consistent with the 
numerical data. This ground-state breaks translational 
symmetry with a unit-cell containing 36 sites, while 36 
sites was the largest sample size that could be studied by 
the exact diagonalizations so far. Therefore, the numer- 
ics might suffer from some serious finite-size effects. In 
particular, the large unit-cell dimer ordering, if it exists, 
will be missed in the current numerical calculations, as 
well as the spinon confinement into magnons. However, 
some gross features of the low-lying spectrum, namely 
the presence of a large number of singlet excitations be- 
low the spin-gap, are not too sensitive to the finite sam- 
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FIG. 10: Monte-Carlo measured spatial spin-spin correlations 
in the disordered phase. The sample size was 20 x 18 x 20 unit- 
cells (each with 6 sites), with h = 0.15, and K$ — Kq — 0.75. 
The type of the site at the origin, and the spatial direction 
are shown in the legend. 



pie sizes used in the numerical studies. (Whether or not 
these singlet excitations have a true gap or are gapless 
in the thermodynamic limit is of course more sensitive to 
finite size). We emphasize that the specific valence bond 
crystal state proposed in this paper is in a confined phase 
where there are no spinon excitations. This is in contrast 
to naive expectations based on the apparent translation 
symmetric ground state obtained in the numerical calcu- 
lations. 

Quite generally, we expect that analysis similar to 
that presented in this paper is applicable to a variety 
of spin models in paramagnetic phases. These include 
the checkerboard lattice, and various other lattices frus- 
trated by the next-nearest neighbor exchange. Then, 
if we were to generalize our hypothesis to other mod- 
els with near neighbor exchange interactions, we would 
expect that the h S> K n limit always describes the ef- 
fective theory for singlet excitations. This theory is a 
quantum dimer model dominated by the kinetic energy. 
A naive guess for the ground-state is that the smallest 
resonating fiippable loops assume the maximum possible 
density, typically leading to a plaquette long range or- 
der. Indeed, such plaquette phases have been observed 
in the Heisenberg model on the checkerboard latticed, 
and the phase diagrams of the quantum dimer model on 
the square^, and honeycomb^ lattices, while a dimer 
RVB phase was found on the triangular lattice^. In this 
paper we demonstrated analytically a plaquette phase in 
the Kagome lattice Heisenberg model. The peculiarity 
of the Kagome lattice is that the maximum density of 
smallest resonating fiippable loops is not sufficient to lift 
the macroscopic degeneracy of the ground-states. This is 
why the singlet gap is extremely small. In other lattices, 
the checkerboard for example, the degeneracy is easily 
lifted, and the singlet gap is fairly significant. 

One of the remaining fundamental questions is whether 
the spin-liquids with gapless singlet excitations (Type II) 
really exist in the quantum spin models, as suggested in 
the Ref. 3,4 . The only known examples which numerically 
exhibit such properties are the Kagome lattice Heisen- 
berg model^, and the triangular lattice antiferromagnet 
with ring-exchangei. In the latter small ring- 

exchange coupling destroys the Neel-order and creates a 
spin-liquid, while at a slightly larger coupling a singlet- 
gap is gradually opened. Armed with our experience on 
the Kagome lattice, we speculate that the case of the 
triangular lattice is also not a Type II spin-liquid with 
gapless singlet states. Probably this phase will not be 
found in other quantum spin models either. 

Another interesting feature that emerged from our cal- 
culations is a rather small (possibly zero) dispersion of 
the excitations in either phase on the dice lattice. The el- 
ementary excitations of both the ordered and disordered 
phases appear dispersionless at least to the sixth order of 
perturbation theory. Although at present we don't know 
whether this persists to all orders, we speculate that it 
might be the case, at least in the disordered phase. In- 
deed, a simple Gaussian Landau- Ginzburg analysis gives 
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a dispersionless spectrum. In the Appendix 151 we study 
the O(N) generalization of the Ising model on the fully 
frustrated dice lattice, and recover a dispersionless spec- 
trum in the large- N limit. In addition, we show that in 
this approximation, the model is always in the paramag- 
netic phase. It is clear that both the Landau-Ginzburg 
approach and the large-A^ approximation have funda- 
mental limitations in this problem. As we have argued, 
the correct phase diagram includes an ordered phase, al- 
beit of a somewhat unconventional kind. 



VII. CONCLUSION 

In this paper we have studied various aspects of the 
fully frustrated transverse field Ising model on a dice 
lattice. This model is useful as a description of the 
physics of singlet excitations below the spin gap of the 
quantum Kagome Heisenberg antiferromagnet. We ar- 
gued that this transverse field Ising model supports at 
least two quantum phases: a plaquette phase where the 
ground state breaks both global spin flip and transla- 
tional symmetries, and a disordered phase with short- 
ranged (Ising) spin correlations. In terms of the origi- 
nal Kagome Heisenberg magnet, the former describes a 
state with a honeycomb pattern of resonating benzene- 
like arrangements of singlet bonds on the hexagonal pla- 
quettes, with a 36-site unit-cell. The elementary singlet 
excitations are gapped with a very small gap and local- 
ized inside the cells of the honeycomb structure. Based 
on comparison with the numerical studies, we suggested 
that this specific large-unit-cell valence bond crystalline 
order describes the ground state of the Kagome Heisen- 
berg quantum magnet. The elementary spin-carrying ex- 
citations are gapped spin-1 magnons. We hope that our 
work stimulates attempts to look for this ordering pat- 
tern in future studies of the Kagome magnet. 

In addition, we argued that due to frustration and un- 
usual geometry the fully frustrated dice lattice transverse 
field Ising model possesses a strong tendency to localize 
its elementary excitations: they are either completely lo- 
calized, or they have an extremely large effective mass 
(narrow energy bands). 
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APPENDIX A: Z 2 GAUGE DESCRIPTION OF 
THE KAGOME QUANTUM MAGNET 

Here we briefly derive a representation of the Kagome 
Heisenberg quantum magnet as a Z2 gauge theory. For 
more in-depth discussion, see referencesii22i2i. Our 



starting point is the Heisenberg model on the Kagome 
lattice: 

H = J^Si- s j ( A1 ) 

(ij) 

where the sum runs over the nearest-neighbor sites. The 
Heisenberg model can be viewed as a special limit of an- 
other model. Consider a theory of hopping spinons f a i 
coupled to the Z2 gauge field which lives on the bonds 
of the Kagome lattice: 

h = -h e °Ti E < h E + h - c -) + 

{ij) (ij) a=U 

+ A l3 (fljl 3 -flfl + h.c.)} (A2) 

The spinon operators obey fermion anticommutation re- 
lations, and the gauge-field operators are Pauli matrices. 
This theory possesses a local Z2 gauge symmetry: if the 
signs of f a i at a particular site i and af t , on all bonds 
emanating from that site are simultaneously changed, 
the Hamiltonian ljA2fl remains invariant. The appear- 
ance of the local symmetry means that the Hilbert phase 
is larger than that of the Heisenberg model, and a con- 
straint per site is needed to project back to the physical 
Hilbert space. This is achieved through a gauge require- 
ment on the physical states: 

G^YlafA-l)^ 1 ^ = -1 (A3) 

The operators Gi are generators of the local gauge trans- 
formations at sites i, and the product is over all bonds 
(ii') emanating from the site i. In the ho — * 00 limit 
of HA2jl . all gauge spins tend to align in the x-direction, 
so that the gauge requirement l|A3(l fixes the number 
of spinons to be one-per-site. The spinon gauge theory 
(|A2I) then reduces to the Heisenberg model at the second 
order of perturbation theory, with J ~ tfj/h , A^/h . 

The next step assumes that the spin-carrying excita- 
tions in the original Heisenberg model are gapped. We 
rely on the experimental and numerical studies which 
show such a gap in the spectrum of the Kagome magnet. 
Thus, the spinon degrees of freedom can be "integrated- 
out" in l|A2|l . leaving behind a pure Z 2 effective gauge 
theory with local terms (we will come back to this step 
later): 

ff ' = -^E4-^EIl4-^EIl4- 

(ij) A A OO 

'v^ElK ••• (A4) 

03 IX 

The terms involving a z operators are the sums over all 
plaquette products. The elementary plaquettes appear at 
the lowest order, but in principle, all connected clusters 
of bond-loops can appear at higher orders. This form 
of the effective Hamiltonian is required by the Z2 gauge 
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symmetry: the sign change of crj/ on all bonds emanating 
from any site i must leave the Hamiltonian invariant. The 
spin-gap results with a local and perturbative nature of 
(|A4|) : the coupling constants K n of the n-bond clusters 
are of the order of ex n , where e is some energy scale, and 
The spinon "integration" also modifies the gauge 
requirement (| A3|) to: 

c i =n°«'=- 1 ( as ) 

i' £z 

The expressions (IA4() and (|A5f) define a pure Z 2 gauge 
theory on the Kagome lattice. We now express this the- 
ory in its dual form. Let vf' v,z be the Pauli operators 
defined on the sites of the lattice dual to Kagome, the 
dice lattice. Their relation to the Z 2 gauge field of the 
Kagome lattice is the following: 

0% = ei m vfv z m (A6) 

n 4 = «f ( a? ) 

plaq. 

In l|A6(l the dual dice bond < Ira > intersects the Kagome 
bond <ij>, while in l|A7(l the dual dice site I sits inside 
the Kagome elementary plaquette appearing in the prod- 
uct on the left-hand side. The numbers e/ m are fixed to 
+1 and —1 in such a way that on every dice elementary 
plaquette the condition of full frustration holds: 

n e '™=- 1 ( A§ ) 



This is the consequence of the gauge requirement (|A5|I . 
One choice of sign-arrangement for e/ m is depicted in 
Fig. n The pure Z 2 gauge theory 1A4II can now be 
rewritten as the transverse field quantum Ising model on 
the fully frustrated dice lattice: 

(hn) I3 1*6 

-K 3+3 J2 (A9) 

The "kinetic energy" terms contain in general the prod- 
ucts of vf on various clusters of the dice lattice sites. 

The Hamiltonian l|A9|l captures only the form of the 
effective theory below the spin gap. The coupling con- 
stants are not known, because the path-integration of 
spinons in (|A2|) cannot be carried out analytically. How- 
ever, we can at least attempt to estimate their rela- 
tive values. For these purposes, we express the original 
Heisenberg model in the slave-particle path-integral, as 
it was done in the Reference^. The action involves the 
spinon Grassmann field f a i living on the sites i of the 
(2+l)D lattice, and the gauge Ising-like field cry living 
on the lattice bonds: 

s = - E & (laifaj + c - c - ] + (A10) 



where the Berry's phase Sb, realizing the projection to 
the physical Hilbert space, is given by: 

e- s -= J] o l3 (All) 

The action l|AT0|) is an exact rewriting of the Heisenberg 
model in the limit £y , Ay <C 1—. We obtain it upon in- 
tegrating out ay, and J At ~ (t?-,A?-), where At — > 
is the imaginary-time lattice spacing used in the path- 
integral (of the Heisenberg model). We can relate the 
action coupling constants iy and Ay to the Hamiltonian 
(IA2|) coupling constants, iy and Ay. In order to sim- 
plify notation, from now on we will use one symbol t to 
represent all the couplings in the action, and t for all the 
couplings in the Hamiltonian. The connection is: 

~ f 1 , on the temporal links (ij) 
13 \ td , on the spatial links (ij) 

where S = \J At /ho — * 0. Now we integrate out the 
spinon field, and obtain a pure Z 2 gauge theory with the 
action: 

S = - log det Aia.j ; iy ) + S B 

= - ^KnUvv + Sb (A12) 
□ □ 

The matrix A above is the matrix which couples the 
spinons in i|A10fl . In the expanded form, this action 
must contain various gauge-invariant products of cry on 
the closed loops, with the loop-dependent coupling con- 
stants K n . The very assumption that this expansion is 
possible and convergent is an essential one, and it re- 
quires a large spin-gap (of the order of J), determined 
at a microscopic scale. A large spin-gap guaranties that 
the effective theory will have only local terms, and sig- 
nificance of various closed loops in (|A12|) will rapidly 
decrease with increasing loop size, making the expansion 
convergent. At this point we are not providing any formal 
estimate of the spin-gap. We simply use some external 
sources to obtain information about it, for example the 
numerics. 

Using the fact that a ay factor always comes together 
with a factor in (|A10fl . we can make an estimate for 
the coupling constants K D (to the lowest order): 

K a - (A13) 

□ 

and this implies that K D can be treated as small num- 
bers, scaling as S n with n being the number of spatial 
links in the loop. In order to formally derive the connec- 
tion between this action and the Hamiltonians (|A4|) and 
(|A9I) . we neglect at this point all the higher order loop 
terms (enclosing more than one elementary plaquette). 
If the action contains only the elementary plaquettes, we 
can readily construct the dual theory, as outlined inil. It 
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takes the form of the Ising model on the fully frustrated 
dual lattice: 



Ira^-lrn^l 



(A14) 



(Ir, 



The Ising gauge-field e; m living on the links of the dual 
lattice is frozen such that its product over any dual el- 
ementary plaquette is —1. The only fluctuating field is 
the vison Ising field Vi- The duality transformation es- 
tablishes connection between the coupling constants K D 
defined for a plaquette and Ki rn defined for the dual bond 
piercing that plaquette: 



tanh Kr 



-2K h 



tanh Ki r , 



-1K n 



(A15) 



For simplicity, assume that there is only one kind of el- 
ementary plaquettes. Then, the Hamiltonian describing 
the dual theory has the form: 



(A16) 



(It. 



The connection between the action (|A14(1 an d Hamilto- 
nian above is: 

KAt = e~ 2Klm , on the temporal dual links (Im) 
hAr = Ki m i on the spatial dual links (Im) 

Now we can estimate the values of the Hamiltonian cou- 
pling constants. The temporal elementary plaquettes al- 
ways have 4 bonds (two spatial and two temporal) , while 
we assume that the spatial elementary plaquettes have n 
bonds. It follows that: 



the path-integral with action i|A10|) to reduce to the 
Heisenberg model. Consequently, h ^> K (for arbitrary 
n > 0). 

This simple analysis suggests that the large-/i limit of 
(IA9I) is a good effective theory for the singlet physics 
below the spin-gap. It also seems natural to expect from 
(|A19|I that the higher-order loop terms, present in the 
more accurate effective Hamiltonian, have the coupling 
constants K n which decay with the loop length n as K n oc 
x n , where x — V J At <C 1. 

And finally, we can ask whether this approach is use- 
ful for extensions of the Heisenberg model. The next- 
nearest and further neighbor exchange interactions are 
reflected in the action IjAlOjl as additional spinon hop- 
ping and pairing terms mediated by the Z2 gauge field. 
The spinons can still be integrated out (provided that 
there is a large spin-gap), and the pure Z2 gauge the- 
ory l|A12l) would effectively live on the more compli- 
cated lattice which has bonds between all pairs of sites 
that host an exchange coupling. The dual theory would, 
in general, be more complicated than the frustrated Ising 
model (the effective lattice may not even be planar) , but 
we could still write an effective Hamiltonian like (|A4|) 
together with the gauge condition (|A5(1 . and expect a 
large-/i limit to be realized. Similarly, the ring exchange 
interactions can be implemented in the action (|A10(I by 
explicitly adding the loop products of the gauge field to 
it. This would yield larger values for the loop couplings 
K n in the pure Z2 gauge theory, possibly making them 
comparable with h. Having in mind the results of this pa- 
per in the small- h limit, this hints that the multiple-spin 
exchange is a good candidate to yield truly fractionalized 
spin liquids. 



J( g—2AT/ m (dual temporal link) 

= tanh K D (spatial plaq.) 
ho 



(A17) 



hAr = if; m (dual spatial link) = 

= — — log tanh K a (temporal plaq.) (A18) 



log[VAr//i ] 

Taking the ratio of i|A17(l and l|A18(l . and noting that 
the Heisenberg model exchange coupling is J ~ t 2 /ho, 
we obtain: 



K 



,/2 



2(JAt) 
|log(JAr)| 



(A19) 



The energy scale of J is also a measure of the spin-gap. 
Since the spinons have to be integrated out, the starting 
path-integral must accurately represent the energy scales 
above the spin-gap. The energy cut-off must be much 
larger than the spin-gap, and hence J At -C 1. This 
is also compatible with request tij , Ay <C 1 , needed for 



APPENDIX B: O(N) LARGE- iV LIMIT 

A relatively simple way to gain some first intuition for 
the quantum Ising model on the frustrated dice lattice 
(J2J) is to generalize it to a limit where the semiclassical ap- 
proximations become exact. Consider the O(N) general- 
ization of spins in the large- N limit. We will demonstrate 
that this Gaussian theory has a dispersionless spectrum, 
and supports only a disordered phase for all values of the 
coupling constants. 

We write the imaginary time path-integral and impose 
the spin-magnitude constraint using a field A of Lagrange 
multipliers. The action obtained from J3J) by keeping 
only the elementary transverse field terms is: 



J dr^—h^2 £imSi T S mT ^ ^ Si 3T 



d 2 S hT 
8t 2 



K 



d 2 S l6T 



2 E^^-E^^- 1 )} ( B1 ) 



In the saddle-point approximation, the Lagrange multi- 
plier field A becomes "homogeneous" , with only two in- 
dependent values A3 and A6 corresponding to the 3 and 
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FIG. 11: A choice of the unit cell on the frustrated dice lattice. 
The site labeling corresponds to the internal indices of the 
matrix in Table ITTTl 



6-coordinated sites. This is the consequence of the Z2 
gauge-invariance, namely the arbitrariness of the choice 
of eim, subject to the condition (J2J). If e; m are chosen to 
be negative on the thick bonds in the Fig.QJb), then the 
action can be diagonalized on the six by six matrices in 
the frequency-momentum space, and the spin degrees of 
freedom can be easily integrated out: 



-iA)+JVtr(iA) (B2) 



Here we used the elementary cell in Fig.^JJ Af is the num- 
ber of space-time sites, A = diag(Ag, Ag, A3, A3, A3, A3), 



and the matrix H qxq w is given in the Table ITTTl The ma- 



trix under the determinant in (IB2|) has three two-fold 
degenerate eigenvalues that have no dependence on the 
momentum. If we relabel iA3 = —m\ and i\§ — — m\, 
they are: 



K; 



V3 



3 2 1 2 

lo + m 3 



u l,2 = -X 



-oj + rrio + m fi ± 



6h 2 + ( 



K 3 - K 6 



-w 1 + mi — m; 



The action then becomes: 



S\ = NJ\f 



2tt 



log(wiw 2 W3) - (2m 6 + 4m 3 ) 



(B4) 



(B5) 



In order for this action to correspond to a physical her- 
mitian Hamiltonian, both m| and m\ must be real num- 
bers. In addition, the path-integral converges only when 
all eigenvalues Vi are real and positive for all frequencies. 
These requirements are satisfied if: 



mi > 



> 



Am\m\ > 6h 2 



(B6) 



The failure to satisfy these conditions by the saddle-point 
values of 771,3 and me would signify an instability. 

Now we vary m\ and m\ to find the minimal action, 
integrate out the frequency, and obtain the following 



saddle-point equations: 



y / 2K 3 m 2 3 
1 



+ K 3 K 6 
-K 3 K 6 



(K 3 - K 6 ) + 2(m '-"^ } 



2(B7) 



where: 



«1,2 




,2\ 2 



K 6 



6h 2 
K 3 K 6 



(B8) 



The conditions (|B6|I directly translate to ai j2 > 0. The 
saddle-point equations are too complicated to be solved 
exactly. If the solution that meets the conditions l|B6|) 
exists for every finite non-zero K 3l Kq, and h, then the 
system is always in the paramagnetic phase. In order to 
prove that such solutions always exist, we can solve the 
inverse problem: assume that m\ and m\ are known, real 
and positive, fix arbitrary values for K 3 and K 6 , and find 
the value of h that satisfies the saddle-point equations. 

A couple of straight-forward algebraic manipulations 
reduce the saddle-point equations (|B7J| to: 



2ml 



K 3 a + 2m 2 - 



= 2- 



1 



2ml \ 2 



/ lm\\ 



= 8AT|a 2 



'm% 



K 3 ) \K 3 

and a — y/ot\ot2- It follows that: 



m| 



K 3 



x 
~2 

2\2 



K, 



(a + x 2 ) 
8K 2 a 2 



(B9) 



(BIO) 



(Bll) 



■a[l- AK 6 {2K 3 + K 6 )a]x + 2K e a 2 = (B12) 



Treating (|B12|) as a quadratic equation for a, we find 
the real and positive solution: 



l + y/l + 8K 6 [2{2K 3 + K 6 )x - 1] 
X AK - 6[2{2K 3 + K 6 )x - 1] 
1 

lor x > 



2(2K 3 + K 6 



(B13) 



We can now study the dependence of (3 = 6h 2 / K 3 K e on 
x. From 1B8I) we obtain: 







K 3 K 6 



(B14) 



and using (|B10|) . (jBTTJ, and l|B13p we express /? as 
a function of x. One can easily check that /3 = for 
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H, 



<3x<?h u 



/ 

/ 









_h„-iq y 

2 C 




2 c 











- 2 W 


2" 


|(1 + e~* 9 « 


\ hi 

) |(- 


-1 + e~' ?y ) 


2 c 






|(l + e iq y) 


h 
2 


K 3 , ,2 
f w 

















h iq y 
2 C 


|(l + e i»v) 





K 3, ,2 
2- W 














2 C 


|(_l + e ^) 










2 







u 


e iqx (l-e iq y) 


_h„iq y 
2 e 













-#- 2 


/ 



(B3) 



TABLE III: The matrix describing the spin-spin couplings in IBljl 



mj = 1/8K 3 , ml = 1/8K 6 , and j3 -> oo for m§ = /iV3/2, 
m| = /i\/3. In between these limiting cases, (3(x) is 
a continuous, monotonously growing function in the re- 
gion x > I/-K3, for all values of K3 and K§. Due to a 
very complicated dependence of (3 on x, we verify this 
by numerical plotting for a variety of K3 and iig val- 
ues, including all important limits. The fact that f3(x) 
is a monotonous function means that for every value of 



h there is a unique solution for x, i. e. for the saddle- 
point equations. The "squared masses" m\ and m\ are 
real and positive for all h, K 3l and K 6 , and vary continu- 
ously with the coupling constants. Therefore, this theory 
does not support a phase transition, and always remains 
in the paramagnetic phase. 



1 C. Zeng, V. Elser; Phys. Rev. B 51, 8318 (1995) 

2 Ch. Waldtmann et al; Eur. Phys. J. B 2, 501 (1998) 



C. Lhuillier, G. Misguich; cond-mat/0109146 (2001) 
P. Sindzingre, C. Lhuillier, J. B. Fouet; |cond-mat/6ll 0283 
(2001) 

M. Mambrini, F. Mila; cond-mat/0003080 (2000) 
A. V. Syromyatnikov, S. V. Maleyev; cond-mat/0304570 
(2003) 

W. LiMing, G. Misguich, P. Sindzingre, C. Lhuillier; 
Phys. Rev. B 62, 6372 (2000) 

Y. J. Uemura et al.; Phys. Rev. Lett. 73, 3306 (1994) 
A. P. Ramirez, B. Hessen, M. Winklemann; Phys. 
Rev. Lett. 84, 2957 (2000) 

A. S. Wills, et al.; Phys. Rev. B 62, R9264 (2000) 
J. B. Marston, C. Zeng; J. Appl. Phys. 69, 5962 (1991) 
S. Sachdev; Phys. Rev. B 45, 12377 (1992) 
G. Misguich, D. Serban, V. Pasquier; Phys. Rev. Lett. 89, 
137202 (2002) 

14 G. Misguich, D. Serban, V. Pasquier; |cond-m at/0302152 
(2003) 

15 R. Moessner, S. L. Sondhi; Phys. Rev. B 65, 024504 (2001) 

16 R. A. Jalabert, S. Sachdev; Phys. Rev. B 44, 686 (1991) 

17 T. Senthil, M. P. A. Fisher; Phys. Rev. B 62, 7850 (2000) 

18 R. Moessner, S. L. Sondhi; Phys. Rev. B 63, 224401 (2001) 



19 R. Moessner, S. L. Sondhi, P. Chandra; Phys. Rev. Lett. 
84, 4457 (2000) 

20 D. Blankschtein, M. Ma, A. N. Berker; Phys. Rev. B 30, 
1362 (1984) 

21 N. Read, S Sachdev; Phys. Re v. Lett. 66, 1773 (1991) 

22 T. Senthil; |cond-mat/0105i04| (2001) 

23 V. Elser; Phys. Rev. Lett. 62, 2405 (1989) 

24 J. Vidal, R. Mosseri, B. Doucot; Phys. Rev. Lett. 81, 5888 
(1998) 

25 L. Balents, M. P. Fisher, S. M. Girvin; Phys. Rev. B 65, 
224412 (2002) 

26 J. B. Fouet, M. Mambrini, P. Sindzingre, C. Lhuillier; 
Phys. Rev. B 67, 054411 (2003) 

27 P. W. Leung, K. C. Chiu, K. J. Runge; Phys. Rev. B 54, 
12938 (1996) 

28 R. Moessner, S. L. Sondhi, P. Chandra; Phys. Rev. B 64, 
144416 (2001) 

29 R. Moessner, S. L. Sondhi; Phys. Rev. Lett. 86, 1881 
(2001) 

30 J. N. Reimers, A. J. Berlinsky, A. C. Shi; Phys. Rev. B 
43,865 (1991) 

31 T. Senthil, M. P. A. Fisher; Phys. Rev. B 63, 134521 (2001) 



